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Abstract — We study physical layer multicasting in multicell 
networks where each base station, equipped with multiple an- 
tennas, transmits a common message using a single beamformer 
to multiple users in the same cell. We investigate two coordinated 
beamforming designs: the quality-of-service (QoS) beamforming 
and the max-min SINR (signal-to-interference-plus-noise ratio) 
beamforming. The goal of the QoS beamforming is to minimize 
the total power consumption while guaranteeing that received 
SINR at each user is above a predetermined threshold. We 
present a necessary condition for the optimization problem to be 
feasible. Then, based on the decomposition theory, we propose 
a novel decentralized algorithm to implement the coordinated 
beamforming with limited information sharing among different 
base stations. The algorithm is guaranteed to converge and in 
most cases it converges to the optimal solution. The max-min 
SINR (MMS) beamforming is to maximize the minimum received 
SINR among all users under per-base station power constraints. 
We show that the MMS problem and a weighted peak-power 
minimization (WPPM) problem are inverse problems. Based 
on this inversion relationship, we then propose an efficient 
algorithm to solve the MMS problem in an approximate manner. 
Simulation results demonstrate significant advantages of the 
proposed multicast beamforming algorithms over conventional 
multicasting schemes. 

Index Terms — Physical layer multicasting, coordinated beam- 
forming, quality of service (QoS), max-min SINR (MMS), 
semidefinite programming (SDP). 



I. Introduction 

With the rapid development of wireless communication 
technology, various kinds of traditional data service, such 
as media streaming, cell broadcasting and mobile TV, have 
been implemented in wireless networks nowadays. As a result, 
wireless multicasting becomes a central feature of the next 
generation cellular networks. Physical layer multicasting with 
beamforming is a promising solution enabled by exploiting 
channel state information (CSI) at the transmitter over the 
traditional isotropic broadcasting. The problem of multicast 
beamforming for quality-of-service (QoS) guarantee and for 
max-min fairness was firstly considered in (Tl. The similar 
problem of multicasting to multiple cochannel groups was 
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then investigated in O . The core problem of multicast beam- 
forming is essentially NP-hard Q . Some other issues such as 
outage analysis and capacity limits were studied in 131. ll4ll. So 
far multicast beamforming has been included in the UMTS- 
LTE / EMBMS draft for next-generation cellular wireless 
services ISl. 161. 

In conventional wireless systems, signal processing is per- 
formed on a per-cell basis. The intercell interference is treated 
as background noise and minimized by applying a predesigned 
frequency reuse pattern such that the adjacent cells use dif- 
ferent frequency bands. Due to the fast growing demand for 
high-rate wireless multimedia applications, many beyond-3G 
wireless technologies such as 3GPP-LTE and WiMAX have 
relaxed the constraint on the frequency reuse such that the 
total frequency band is available for reuse by all cells in 
the same cluster. However, this will cause the whole system 
limited by the intercell interference. Consequently, cooperative 
signal processing across the different base stations has been 
identified as a key technique to mitigate intercell interference 
in the next-generation wireless systems. 

The goal of this paper is to investigate multicast beam- 
forming for cooperative multicell networks. The base sta- 
tions cooperate with each other and design their transmit 
beamformer s in a coordinated manner. In general, there are 
two cooperation scenarios in multicell networks. In the first 
scenario, different base stations are fully cooperative and act 
as "networked MIMO (multiple-input multiple-output)" (e.g., 
171 . ISl). Namely, they coordinate in the signal level, i.e., data 
information intended for different users in different cells is 
shared among the base stations. Clearly, networked MIMO 
needs tremendous amount of information exchange overhead. 
In contrast, in the second scenario, the base stations are only 
required to coordinate at the beamforming level which needs 
rather small information sharing (e.g., 0). In this paper, we 
focus on the latter case by only allowing beamforming level 
coordination. 

We first formulate the problem of multicell multicast beam- 
forming as sum power minimization subject to the constraint 
that the received signal-to-interference-plus-noise ratio (SINR) 
of each user is above a threshold. This problem is referred to 
as QoS problem. It is known that the QoS problem in single- 
cell scenario is always feasible. However, this is not so in 
multicell networks due to the intercell interference, especially 
when the channel condition is bad or the SINR target is 
very stringent. We first present a necessary condition on the 
optimization problem to be feasible. Then we propose a novel 
distributed algorithm to solve the multicell multicast coop- 
erative beamforming design. More specifically, based on the 
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decomposition theory ifTOl , we introduce a set of interference 
constraint parameters and decompose the original problem 
into several parallel sub-problems. Since the sub-problems are 
nonconvex and NP-hard, each base station then applies the 
semidefinite relaxation and independently solves its own sub- 
problem. The interference parameters are updated based on 
a master problem. Simulation results show that the algorithm 
converges in only several iterations. 

In view of the user fairness issue, we also consider the multi- 
cell multiast beamforming design with the goal of maximizing 
the minimum SINR (MMS) among all users under individual 
power constraints. This is called MMS problem. By linking 
the MMS problem with a weighted peak power minimization 
(WPPM) problem, we propose an efficient algorithm to find 
the near-optimal solution. Previously, the authors in ifTTIl |[T2l 
also studied the multicell multicast beamforming problem with 
the objective to maximize the minimum SINR but subject to a 
total sum power constraint across all base stations. Hence, the 
problem considered in ifTTIl |[T2ll is very similar to the multicast 
beamforming problem in single-cell multi-group systems as in 
llJl. Authors in [13] also considered the similar problem under 
per-base- station power constraint but allowed data sharing 
among base stations, which is similar to the "network MIMO" 
and requires a large amount of information exchange overhead. 
Thus, our beamforming design for max-min fairness is more 
practical as we consider individual peak power constraint 
on each base station and only allow coordination in the 
beamforming level. 

The rest of the paper is organized as follows. In Section II, 
the system model for multicell multicasting is presented. Sec- 
tion III considers the beamformer design in the QoS problem 
and develops a novel decentralized coordinated beamforming 
algorithm. In Section IV, we consider the max-min SINR 
beamforming problem under individual base station power 
constraints. Section V provides simulation results. Concluding 
remarks are made in Section VI. 

Notations: R and C denote the real and complex spaces. 
The identity matrix is denoted as I and the all-one vector is 
denoted as 1. For a square matrix S, S h means that S is 
positive semi-definite. denotes its element in the ith row 

and jth column. {■)^, {■)^, {-Y and Tr{-} stand for transpose, 
Hermitian transpose, Moore Penrose pseudoinverse and the 
trace respectively. \x\ denotes the absolute value of the scalar 
X and llxll denotes the Euclidean norm of the vector x. 



II. System Model 

Consider a multicell multicast network comprising N cells 
and K mobile users per cell as shown in Fig. 1. The base 
station in each cell is equipped with Nt antennas and ev- 
ery mobile user has a single antenna. Let h^^^ denote the 
frequency-flat quasi-static 1 x TV^ complex channel vector from 
the base station in the ith cell to the kth user in the jth cell, and 
denote the NfXl multicast beamforming vector applied to 
the base station in the ith cell. We define the complex scalar 
Si as the multicast information symbol for the users in the ith 
cell. The discrete-time baseband signal received by the kth 




Fig. 1. Multicell multicast network. 



user in the ith cell is given by 



N 



yi,k 



.WjSj 



(1) 



where Zi^k is the additive white circularly symmetric Gaussian 
complex noise with variance (jf^/^ on each of its real and 
imaginary components. In (1), the second term is the intercell 
interference. 

Based on the received signal model in (1), the performance 
of each user can be characterized by the output SINR, defined 
as 



SINR, 



(2) 



Notice that each user in the considered multicell multicast 
network only suffers from the intercell interference, which is 
different from multicell unicast systems where both inter-cell 
interference and intra-cell interference exist. 

In practical scenarios, the channels from a base station to 
different users, which may or may not belong to a same cell, 
can be correlated, especially when these users are in close 
proximity to each other. For the users belonging to the same 
cell, we call the correlation between h^^^ and h^^^, for k ^ 
k' as the intracell-user channel correlation. For the users from 
different cells, we call the correlation between h.f- ^ and h.f-, ^ 
for j ^ j' as the intercell-user channel correlation. 

III. QoS Beamforming 

The coordinated QoS beamforming design is to minimize 
the total energy consumption of the system while maintaining 
a target SINR for all users by properly designing the beam- 
formers at each base station. This is formulated as: 



N 



P(7) : min ^ ||wi 



s.t. 



N 



(3) 

>7i,Vi,fc (4) 



where 7 = [71, 72, Tat]^ is the target SINR vector with 
each element 7^ being the target SINR value to be achieved 
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by the users in the ith cell. Since the base station transmits 
a common information in a multicast manner, the information 
rate for the users within one cell is the same, and hence we 
set a common SINR target for all the users in the same cell. 

A. Feasibility analysis 

Due to the SINR constraints, the QoS problem in ^ is 
not always feasible, which is similar to the multiuser unicast 
scenario {lAj . In order to verify its feasibility, we need to show 
whether there exists a set of beamformers {^/^^}fLl for a given 
7 such that 



min SINR,,fc>7,,ViG{l,2,...,A/-} 

k 



(5) 



For simplicity, the SINR targets for different cells are assumed 
to be the same. Then the above condition boils down to the 
following 

min SINRi,/e > 7 (6) 

i,k ' 

For the kth user in each cell, we combine their channel vectors 
as follows 



1,1, 



."l,Ar,/e 



"2,1,/c 



"2,Ar,/c 



"Ar,2,fe 



k = l,2,...,i^ 



(7) 

where each H/^ is an TV x (TV x TV^) matrix. The following 
lemma provides a necessary condition for the QoS problem to 
be feasible. 

Lemma 1: Given the SINR target vector 7 = 7 • 1, if the 
problem ([3]) is feasible, then the SINR target 7 should satisfy 
the following condition 



7 < min 



rank(Hi) rank(Hx) 
TV-rank(Hi)''" ' TV - rank(HK) , 



(8) 



Proof: Please refer to Appendix A. ■ 
From Lemma 1, we can see that the intercell-user channel 
correlation has a negative impact on its feasibility. More 
specifically, if some of the rows in H/e, say the ith row and i'th 
row for i i\ are correlated, which means that the channels 
between the kth user in the ith cell and the kth user in the 
i'th cell are correlated, then the rank of H^^ could be less than 
TV and as a result, the SINR threshold 7 will have a finite 
upper bound. On the other hand, if the channels of users in 
different cells are independent, we have that the matrix H/e is 
full rank with probability one, then the corresponding upper 
bound is infinite, which means no constraint on 7. Since the 
problem © is NP-hard [2], determining its feasibility is not 
an easy job. The lemma gives a necessary condition from the 
perspective of channel correlation. In the following, we only 
consider the problem ^ when it is feasible. 

B. Decentralized coordinated beamforming 

A desired feature of coordinated beamforming in multicell 
networks is that the base station at each cell can implement 
its beamforming design locally ifTTIl . ifTSl . This is due to the 
constraint in practical systems that the backhaul channel has 



limited capacity. In this subsection, we propose a decentralized 
scheme for implementing the multicell multicast QoS beam- 
forming. It is assumed that each base station in the network 
only has the channel knowledge of the mobile users within its 
own cell. 

The distributed algorithm to problem P(7) can not be easily 
obtained primarily because all the beamformers are coupled 
together in the constraint ©. According to the decomposition 
theory ifTOl and inspired by ifTSll , we introduce a set of slack 
variables Tij^k denoting the constraint of the interference from 
the ith base station to the kth user in the jth cell. Then the 
problem P(7) can be reformulated as 



P(7,r) : 



mm 



S.t. 



N 

El 



N 



< F 



(9) 



>7„Vz,fe (10) 



Here, the real-valued vector T is defined as follows 

r = (ri.2.i,...,ri.2.K,...,r^,^-i.ifreiEi(^(^-i)^)><^ 

(12) 

The introduction of T is similar to the concept of interfer- 
ence temperature (IT) in cognitive radios and hence we refer 
to it as IT vector. It is now observed that the constraints in ([TOb 
are all decoupled. Unlike ifTSl , our problem is still nonconvex 
and NP-hard. We cannot solve it by exploiting its dual problem 
since the strong duality does not hold. In the following we take 
the semidefinite relaxation (SDR) approach, the optimality of 
which shall be discussed later in this section. Introducing new 
variables {W^ = w^wf }^^, the relaxed problem of P(7, F) 
becomes: 



N 



(13) 



s.t. Tr{Hi,i,fcWi} > 7ie^,,fcr + ^^(Tlkyi, k (14) 
Tr{Hi,,- fcWO < e^^- fcF, Vi,i ^i,k (15) 



W. ^0 



(16) 



Here, for notation convenience we have introduced the 
(TV {N — 1) K) X 1 direction vectors eij^k and defined 

'^i,j,k = ^i,j,k^i^j^k' 

For a pre-fixed IT vector F, since the constraints have 
been decoupled, problem Pi(7,r) can be decomposed into 
TV parallel subproblems. The ith subproblem is as follows 

P,^^^(7,r): minTr{Wj (17) 

Wi 

s.t. Tr{H,,,,fcW,} >7,e^,,fer + 7,af,fe,V/c (18) 
Tr(H,,,- W,) < e^^. ^r, Vj ^i,k (19) 
WihO (20) 

It is easily observed that in each subproblem, the base station 
only needs the local channel state information. Specifically, 
the ith cell's base station needs just the channel vectors 
h^^^, Vj = 1 • • • TV, V/c = 1 • • • i^T. Since it is also convex, the 
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optimal solution can be obtained efficiently. Having solved the 
subproblems, we then define the master problem in charge of 
updating the IT vector F 



formulated as below 



>mas 



(7): mm P(r) 



(21) 



N 



where P (T) = Y.Pii'^) with i^*(r) being the optimal 

solution of problem P^^^(7,r) for a given F. This master 
problem can be solved iteratively using a subgradient pro- 
jection method. Denote g G IR(^(^-i)^)xi as the global 
subgradient of P (F) at F. The following theorem suggests 
that g can be obtained from each base station. 

Theorem 1 : The global subgradient g of P (F) in the 
master problem P^^^(7) is given by 



N 



(22) 



where g^ is the subgradient of P/(F). 

Proof: Please refer to Appendix B. ■ 
According to Theorem 1, every base station first solves 
its own subproblem, gets the subgradient vector g^ and then 
broadcasts it to other cells. Upon receiving all the N subgra- 
dient vectors g^'s, the base station in each cell will sum them 
up to get the subgradient g and update the IT vector F as 
below 



F(n + 1) 



F(n) - /i(n) 



1 + 



|g(n) 



(23) 



where n denotes the iteration index and /i is the step size. Here, 
[•]+ denotes the projection onto the nonnegative orthant. Since 
the problem Pi(7,F) is convex, this distributed algorithm is 
guaranteed to converge and converge exactly to the optimal 
solution of Pi(7,F). Just like the steepest descent method, 
the choice of step size affects the convergence properties of the 
iterative algorithm, such as the speed and the accuracy. Here 
we just take the simple diminishing step, i.e., /i(n) = sj ^Jn, 
where 5 > is the initial step size. 

When the algorithm converges, each base station can get 
its own beamformer from its resulting matrix WJ". We now 
discuss how to extract the beamformer vector from each 
W* and its optimality. Since the original problem P(7) 
is NP-hard, generally there is no guarantee that an algo- 
rithm for solving the relaxed SDP problem will give the 
desired rank-one solution. If the WJ^ is rank-one, then the 
base station applies the eigen-value decompsotion (EVD) to 
as = \\^\^\^ and takes = as 

the optimal beamformer. Otherwise, it first generates a set 
of candidate beamforming vectors {w-} by randomization 
method. Specifically, one can perform EVD on WJ" to get 
WJ" = U^E^Uf^ and then generate the /th candidate vector 
w- as w- = U^XlJ^^v/, where ~ CA/'(0, 1). Then based on 
its own subproblem P^^^(7, F), the base station needs to do 
scaling to get the beamforming vector. The scaling problem is 



L^(7, F) : min di 



ai > 



(24) 

(25) 
(26) 
(27) 



The above scaling problem is a linear programming prob- 
lem about the one-dimension variable oci and can be solved 
very fast. Denote the minimum scaling ratio as ol\ and the 
associated vector as w^, then the corresponding beamformer 
is = ^J~otl • wj". Both the randomization and scaling 
procedures are implemented locally by each base station and 
do not break the distributed nature of the algorithm. Notice 
although the proposed algorithm cannot guarantee the optimal 
solution for the original problem P(7), it is seen from the 
simulation results in Section V that this algorithm achieves 
the optimal solution of the original problem in most cases, 
especially when the network size is small. 

Finally, the decentralized algorithm is summarized below. 

Algorithm 1: Decentralized algorithm for QoS beam- 
forming 

• Step 1. Initialize the IT vector F(0) with certain values, 
e.g., F(0) = 1 and set the iteration number n = 0. 

• Step 2. Each base station locally solves its sub-problem 
([TtI) and broadcasts the resulting g^ (via the backhaul 
signaling) to the other — 1 cells' base stations. 

• Step 3. Based on the reception of the subgradient vectors 
from all other cells, each base station updates the IT 
vector F(n) according to (l23l) . 

• Step 4. set n = n + 1 and go to Step 2 until meet the 
stopping condition. 

• Step 5. For each base station i, if is rank-one, 
then the optimal beamformer is the eigenvector of 
W^; Otherwise the base station implements the random- 
ization and scaling procedures to get the near-optimal 
beamformer. 

The main information exchange in this decentralized beam- 
forming scheme is the real- valued subgradient vectors. During 
each iteration, for each base station, it should broadcast its 
subgradient vector g^, which only has N x K nonzero entries. 
The sum signaling overhead among the base stations in one 
iteration is thus O (^N'^K). Denote the number of iteration 
times as A^5, then the total signaling is O (^N^N'^K). Through 
the simulations shown in Section V, we can see that the 
algorithm converges very fast and can achieve major part of the 
beamforming 's gain after only several iterations. Furthermore, 
this distributed algorithm can also work for the conventional 
interference channel, which is a special case (i^T = 1) of the 
multi-cell multicast networks. 

IV. Max-Min SINR Beamforming 

In this section, we study the beamforming problem for the 
MMS scheme under individual base station power constraints. 
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The problem is formulated as follows 
S(p) : max min 



s.t. 



(28) 



(29) 



where p = [pi,P2, •.•,PAr]^ is the power constraint vector for 
base- stations in the network. 



A. Connection with power optimization 

Define a weighted peak power minimization problem as 



Q(7,p): min max — ||w^|p 



s.t. 



N 



(30) 

>'^,Mi,Mk (31) 



where 7 is the common SINR target for all cells. It can be seen 
that the two problems S(p) and Q(7,p) are connected with 
each other through the power vector p. We use the notation 
7 = S(p) to stand for the optimal objective value of problem 
S(p), which means that the maximum worst-case SINR for 
S(p) is 7. For the problem Q(7, p), we denote the associated 
optimum value as p = Q(7,p). Then the following theorem 
tells the relationship between these two problem: 

Theorem 2: The SINR optimization problem of (|28]) and 
the power optimization problem of (l30l) are inverse problems: 



7 = S(Q(7,p)-p) 
l = Q(S(p),p). 



(32) 
(33) 



Proof: We first prove (1321) by contradiction. Let p and 
{^i}iLi be the optimal solution of Q(7,p), and 7 and 
be the optimal solution of S(p • p). We assume that 
7^7. Then if 7 < 7, then we can choose {w^}^^ as the 
solution for S(p • p), which provides a larger objective value 
7. This is a contradiction of the optimality of {w^}^^ for 
S(p • p). Otherwise, if 7 > 7, then we can find a constant 
c < 1 to scale the solution set {w^}^^ while still satisfying 
the SINR constraints of problem Q(7,p). Since {w^}^^ 
satisfy the power constraints in S{p • p) which means that 
maxvi^||wi|p = p, the resulting set {cwi}f^^ achieve a 
smaller objective value (weighted peak base- station power) 
than p, which contradicts the assumption that {w^}^^ is the 
optimal solution of Q(7, p). Thus we must have 7 = 7. 
The proof of ([33l) is similar and therefore omitted. ■ 
In addition, we find numerically that the optimal objective 
values of both the two problems are monotonically non- 
decreasing in the constraint parameters p and 7. Such finding 
is reasonable since with more power on each base station, the 
larger SINR can be achieved, and vice versa. 

A similar theorem on inverse property for single-cell multi- 
casting has been proved in O. Unlike ||2l, our max-min SINR 
problem for multicell multicasting cannot be solved by directly 
solving the corresponding QoS problem since each base station 
is subject to an individual power constraint. Instead, we 



show the this max-min- SINR problem with multiple power 
constraints can be solved by solving another weighted peak 
power minimization problem. 

B. Inversion-property-based algorithm 

Since the problem S(p) is non-convex, similar to section 
III, we apply the semidefinite relaxation and get the relaxed 
problem Si(p) as below 

Tr{H,,,,fcW,} 



Si(p) : max min — — 



(34) 



ETr{H,,i.feW,} + a2, 

S.t. Tr{W,}<p„Vz (35) 
WihO.Wi (36) 

The non-convex rank-one constraint has been dropped. We first 
get the optimal solution of Si(p) based on Theorem 2. The 
inverse problem of problem Si(p) is just the relaxed version 
of Q(7,p), which is defined as follow 



Qi(77P) • ™n iiiax — TrjWil 



(37) 



N 



S.t. Tr{H,,,,feWj > 7^Tr{H,-,,fcW,} + 7a^^^,Vz,Vfe (38) 

WihO.Wi (39) 

Introducing a slack variable x, we rewrite this problem in a 
more elegant way 



Qi(7,p) : min x 



(40) 



N 



S.t. Tr{H,,,,fcW,} > 7^Tr{H,-,,fcW,} + 7a^,^,Vi,Vfe (41) 



—Tr{Wi} <x,\fi 
Pi 



(42) 
(43) 



Then we can solve Si(p) by iteratively solving its inverse 
problem Qi(7,p) for different 7's. Notice that Qi(7,p) is 
a SDP problem with strong duality and thus can be solved 
efficiently using the interior method. Due to the inversion 
property, if 7q is the optimal value for Si(p), then the optimal 
value for the problem Qi (7q , p) should be equal to 1. With the 
non-decreasing monotonicity, we can find the optimal value 7q 
efficiently by a one-dimensional bisection search over 7. When 
we get the optimal solution of Si(p), based on the resulting 
matrices {W^}f^, we apply the EVD or the randomization 
and scaling to obtain the final beamformers for all the base 
stations. Finally, the beamforming algorithm is summarized 
below 

Algorithm 2: Inversion-property-based algorithm for 
MMS beamforming 

• Step 1. Initialize the interval [L^U] which contains 
the optimal value 7q of Si(p), e.g., L = O^U = 
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maxvi,v/c Pi • ||hi,i,fe|P/c^^, and set the iteration number 
n = 0.' 

• Step 2. Set 7 = {L -\- U)/2, and solve the problem 
Qi(7,P)- 

• Step 3. If the optimal value Pq of Qi (7, p) is larger than 
1, set /7 = 7; Otherwise, set L = 7. 

• Step 4. Set n = n + 1 and go back to Step 2 until meet 
the stopping condition. 

• Step 5. If WJ" is rank-1 for all i, we can obtain the 
optimal solution for S(p) by EVD; Otherwise, the central 
controller does randomization and scaling to obtain the 
approximate solution. 

Note that, in each iteration, we need to solve the weighted 
peak power minimization problem, which is the main differ- 
ence from the algorithm in O . Although this algorithm gives 
an approximate solution for S(p), from the simulation results 
in Section V, we can see that it obtains the optimal solution 
in most cases. 

V. Simulation Results 

In this section, we provide numerical examples to illustrate 
the performance of the proposed multicell multicast beam- 
forming designs. Within each cell, the channel is assumed as 
the normalized Rayleigh fading channel, i.e., the elements of 
each channel vector are independent and identically distributed 
circularly symmetric zero-mean complex Gaussian random 
variables with unit variance. For the intercell channels, i.e., 
the channels from the ith cell's base station to the users in the 
jth cell (j ^ i), we introduce the average large-scale fading 
ratio £, < £ < 1. A big e means large intercell interference, 
which often occurs when the users are at the boundary of 
each cell. Here, we set £ = |. A common noise variance is 
set to be cr|^ = 1, Vi, /c. Throughout this section, we use the 
notation 'W — K — Nt' to describe the system configuration, 
which means that the network has N cells with K mobile 
users in each cell and each base station has Nt antennas. For 
all simulations, 200 channel realizations are simulated and 100 
Gaussian randomizations are generated if the randomization 
method is needed. 

A. Performance comparison with existing schemes 

In this subsection, we consider the performance of the 
proposed algorithms for the two transmission schemes. For 
the QoS scheme, we assume that the SINR targets for the 
users in different cells are the same for simplicity. We first 
illustrate the convergence of the algorithm. Fig. 2 plots the 
power consumption during each iteration at different target 
SINR for the (2 — 2 — 4) network. The initial step size is /i = 1 
and the diminishing step size is /i(n) = ii/^/Jn). Since the 
problem Pi (7) is convex, the proposed distributed algorithm 
always converges to its optimal value. The simulation result 
validates it. Also we can see that at the first few iterations, the 
algorithm converges very fast and achieves the major part of 
the optimal value. 

In Fig. 3, we illustrate the efficiency of the SDR approach 
by comparing it with the lower bound, which is the solution 




40 
20 
0, 



>Q o o o o o o o o o o 4 > 

SINR target = 6:dB 



SINR target = 6:dB 

o o <gj o O Cp o o o o 



5 10 
iteration number 



15 



Fig. 2. Convergence behavior of the decentrahzed algorithm for the QoS 
scheme. 
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Fig. 3. Comparison of the proposed algorithm with its lower bound for the 
QoS scheme. 

of problem Pi (7) in ([T3l) where the rank-one constraint is 
relaxed. It can be seen that although we take the semidefinite 
relaxation, the proposed coordinated beamforming algorithm 
obtains the optimal performance in most cases. 

We now compare the performance of the proposed SDR- 
based QoS beamforming in Algorithm 1 with two conven- 
tional multicell beamforming algorithms: multicell block di- 
agonalization (M-BD) 1 19] and layered signal-to-leakage-plus- 
noise ratio (L-SLNR) |20|. Though these two algorithms were 
proposed for multicell networks with unicast traffic, they can 
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Fig. 4. Comparison on feasibility of the different multicast beamforming 
algorithms for the QoS scheme. 
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be easily extended to multicasting. More specifically, for the 
M-BD algorithm, each base station chooses the beamforming 
vector which lies in the null space of the channels from the 
other — 1 base stations, i.e.. 



Wi e Null 



(«') 



(44) 



where H, 



hi+i,i,i, hi+i,i,K, hAr,i,i, hN,i,K] and Null(.) 
stands for the null space of a matrix. Then the QoS problem 
reduces to a power allocation problem and we can solve it to 
obtain the resulting power for each base station. For the L- 
SLNR algorithm, each base station chooses its beamforming 
vector as 



Wi oc max.eigenvector 



N K 



-^ -1 



K 



Eh. 

.1=1 



(45) 



which means that should have the same direction as the 
eigenvector corresponding to the largest eigenvalue of the 
above matrix. We can similarly get the resulting beamforming 
vectors for the base stations by solving the reduced power 
allocation problem. As a performance benchmark, open-loop 
space time block coding (STBC) without requiring any CSI 
at each base station is also simulated. Fig. 4 shows the 
comparison on feasibility 0. Here, the feasibility percentage 
is obtained by counting the number of trials among the 200 
channel realizations that each algorithm is able to find the 
solution to meet all the constraints. It can be seen that both 
the open-loop STBC and the L-SLNR beamformer almost 
do not work when SINR target 7 is large. This is expected 
as the open-loop STBC does not make any use of channel 
state information and serves purely as isotropic broadcasting. 
The L-SLNR beamformer, on the other hand, only tries to 
maximize the SLNR from the transmitter perspective and 
cannot guarantee the SINR maximization at the receiver end. 
From Fig. 4 it is also seen that the M-BD beamformer and the 
proposed coordinated beamformer are always feasible in the 
considered SINR target region. Here, the reason that the M- 
BD beamformer can work well is that the considered network 
is an interference-limited one and M-BD can null out all 
the interference for each user. The good performance of the 
proposed algorithm is expected as it can obtain the optimal 
solution in most cases under this particular setting (see Fig. 
2). Therefore, in Fig. 5 we only compare the energy efficiency 
of the proposed algorithm with M-BD algorithm. We can see 
that the M-BD algorithm consumes much more power than 
the proposed algorithm. In particular, at target SINR of lOdB, 
M-BD consumes 3 dB more power in the (2 — 2 — 4) system 
and 4dB more power in the (3 — 2 — 6) system. Notice that 
the M-BD algorithm requires that the number of transmitting 
antennas at each base station should be larger than the total 

^Although we have assumed that the beamforming problem is feasible, each 
specific algorithm may not be able to find the beamformers to meet the SINR 
target due to its sub-optimality. 
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Fig. 5. Comparison on energy efficiency of the different multicell multicast 
beamforming algorithms for the QoS scheme. 
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Fig. 6. Comparison on the minimum SINR of the proposed algorithm with 
its upper bound for the MMS scheme. 

number of receive antennas at all users. However, our proposed 
algorithm does not have this requirement. 

For the max-min SINR beamforming scheme, we assume 
the power constraint for each base station is the same. We first 
compare the the performance of the proposed Algorithm 2 
with its upper bound in Fig. 6. Similar to the QoS scheme, we 
can see that the gap from the upper bound is very small and it 
achieves the optimal solution in most cases. Fig. 7 compares its 
performance with the existing algorithms. We can see that our 
proposed algorithm significantly outperforms the other ones 
over a large range of individual power constraint parameter. 
In particular, at the individual power constraint of lOdB in 
the (3 — 2 — 5) system, the worst-case SINR achieved by the 
proposed algorithm is 6dB higher than L-SLNR, 8dB higher 
than M-BD and 9dB higher than open-loop STBC. From Fig. 
7, it is also found that the M-BD algorithm performs the worst 
when the per base station power is small but is superior to the 
L-SLNR algorithm at large per base station power. 

B. Ejfects of channel correlation 

In practical scenarios, the channels among different users 
may be correlated with each other, especially for the users 
who are very near to each other in the geographical position. 
In order to model the correlated channels, we use the following 
Kronecker model (Ell ES- 

-Rkro = C^/^G, (46) 
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Fig. 8. Feasibility of the proposed algorithm with correlated intercell-user 
channels for the QoS scheme. 



where tLkro is a channel matrix whose rows are correlated 
with each other and G is a matrix with i.i.d. circularly 
symmetric Gaussian entries with zero mean and unit variance 
(This is for intracell channels by default. If it is intercell 
channel, the variance is e^). We model channel correlation 
matrix C as a Hermitian Toeplitz matrix with exponential 



entries [C]i^- 



. Here, r can be seen as the 



correlation ratio and < r < 1. 

We first investigate the effects of intercell-user channel 
correlation on the feasibility of the QoS problem in Fig. 8. 
The correlation ratio r is set to be 0.5, 0.7 and 0.9, where 
r = stands for the independent channel. We can see that 
when the intercell-users' channels are correlated, the feasibility 
decreases, which justifies our statement in Section III- A, i.e., 
the intercell-user channel correlation has a negative impact on 
the feasibility of the QoS problem. 

Fig. 9 shows the effects of intracell-user channel correlation 
for the QoS beamforming scheme. The correlation ratio r is 
also set to be 0.5, 0.7 and 0.9. From the results, we can see that 
when the intracell-user channels are correlated, it is helpful for 
the system. More specifically, the consumed power becomes 
less. We have observed the similar results on the usefulness of 
intracell-user channel correlation on the MMS scheme, which 
are ignored here due to page limit. 



VI. Conclusion 

This paper considered two coordinated mulicast beamform- 
ing schemes for multicell networks. For the QoS scheme, we 
provided a necessary condition for the beamforming problem 
to be feasible when the network shares a common SINR 
target and proposed a decentralized algorithm to implement 
the coordinated beamforming. For the max-min SINR scheme, 
we considered individual base station power constraints and 
also proposed an efficient beamforming algorithm. Besides, we 
also investigated the impacts of intercell-user and intracell-user 
channel correlation on the multicast network. 

There are several other issues to be further investigated for 
multicell multicast beamforming. First, finding the sufficient 
condition for the feasibility of the QoS problem remains 
open. Second, a decentralized algorithm to implement the 
beamforming design for the MMS scheme is to be designed. 
Last but not least, it is also interesting to consider robust 
beamforming when the channel state information at each base 
station is not perfect. 



Appendix A 
Proof of Lemma 1 

We first define the beamformer matrix W which is full-rank 
as below 



W 



wi 

W2 



(47) 



••• WN_ 
Based on the fact that SINR^^/c < SIR^^/c, we have 

SINR, , 



N 



<J 

— N 

Y. 



- 1 

Vi,k 



(48) 
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where we define the r]i^k as below by 



|h?, 



|2 



, = ^M.^^^l L{Wu r, A) = Tr{(W,)} 



^ n 



Tr{H,,,,fcWJ - e"T 



(56) 



the second equation in m is because © and (StJ. We also ^j^^^ ^j^^ ^^^j ^^^^^j^^ ^^^^^ 
know the monotonicity of f{x) = 1/ [l/x — 1) when x < 1. 

Comparing it with (|48l) , we can bound the SINR^ ^ through t /t^ a n _ • j /^^r a n 
1 • 1 • T 1 • 1 1 i' • V 5 ^/ — min W ^ , -I- 7 ^ ) 

bounding the associate rji^k- Let the singular value decomposi- 

tion(SVD) of H/e W be denoted as H^W = UfeXIfe V^, where / k n k 

both Ufe and Vfc are K x Vk quasi-unitary matrix, i.e., their = ^ h,i,k^J^i^k ^i.j.k^J.j.k 1 ^ + /i(A) (57) 

columns are orthogonal with each other. 5]^ is an r/c x r/c \ /c=i ji-i k=i 

diagonal matrix whose elements are the singular values and 
Tk = rank(H/eW). Then where 



|2 



^.,fe = ^^%^f^, ^ = 12.--7V,fe = 12---i^ (50) /^W 



where Uk,i and ^rk,^ are the /th columns of \Jk and V^. ( Tr{W^} ^ 

According to Cauchy-Schwarz inequality, we have 



i,k 

N K 



-Tr{H,,,,fcW,}-af^ 

7i 



|<,5],v.,|' < |K,5],f • ||v,,f A,,,,Tr{H,,,,Wj (58) 



K^E^.u,,,) K^^v,,,) (51) 



jVi k=i 



Since vf .v/e,i = [(H/eW)^(H/eW)]i,i, we conclude that Since the problem Pf^^(7,r) is convex, then the strong 

duality holds which means that 

7/,,fe < [(HfcW)t(HfeW)],,, (52) 

i^*(r)=maxd,(r,A) (59) 
Thus we have a^o 

N N 

^ 1 ^ 1 v^rz-r-r ,,rx + /-rT ^^T\^ Dcnotc A"^ as the optimal Lagrange multiplier for the dual 

. n,K - n,K - K J \ K problem, then we have 

= lTr{(H.W)tH.W} = < !^ ^*(r) = d.(r, V) 

(53) I K N K \ 

= E K^,Ai,k r + /.(A*) (60) 

Further we can get .^^^^^ j 

inm = mm {m.in 7?^,^ } Define & as 

" ' ^ ' ' ' ' g.^E^U.e...-EE%^-...^' (61) 

So if the problem © is feasible, then the minimum SINR k=i j^i /e=i 

should be larger than the threshold 7. Plugging (l54l) into (l48l) , 

we can conclude that we have 

/ • cTXTo rank(Hi) rank(HK) 1 d^^t-^ «.^t-i /-^x^^ 

hk I A/^ — rank(Hi) A/^ — rankiH^) I rr , rr- , rr, ^ 

^ ^ ^ ^ (Vs) =gf(r-r) + gfr + /,(A*)<gf(r-r) + />*(r)(62) 



which completes the proof of Lemma 1 . 

Appendix B 
Proof of Theorem 1 



which is equivalent to that 

P:{f)>P*{T)+gf{f-T). (63) 



We begin by computing the subgradient from each It means that is the subgradient of P/(r) and obtained for 
subproblem. The Lagrangian of the ith subproblem is given the ith subproblem. 
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In the same way, we can compute the global subgradient g 
of P"^^^(7,r) as below 

N K N N K 

i—1 k—1 i—1 j^i k=l 



N 

E 



K 



i=l \/c=l 
N 



N K 



(64) 



which completes the proof of Theorem 1 . 
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